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Abstract 



In this paper, we provide algorithms to rank and unrank certain degree-restricted classes 
of Cayley trees. Specifically, we consider classes of trees that have a given set of leaves or that 
have a fixed number k of leaves. Using properties of a bijection due to Egecioglu and Remmel 
[3] , we reduce the problem of ranking and unranking these classes of degree-restricted trees 
to corresponding problems of ranking and unranking certain classes of set partitions. For 
fixed k, the number of Cayley trees with n vertices and k leaves grows roughly as n! and 
hence the ranks have 0(nZog 2 ("-)) bits. Our ranking and unranking algorithms require at 
most 0(n 2 ) comparisons of numbers y < n plus 0{n) operations of multiplication, division, 
addition, substraction and comparision on numbers x of length 0(nlog(n)). 

1 Introduction 

In computational combinatorics, it is important to be able to efficiently rank, unrank, and 
randomly generate (uniformly) basic classes of combinatorial objects. A ranking algorithm for a 
finite set S is a bijection from S to the set {0, • • • , — 1}- An unranking algorithm is the inverse 
of a ranking algorithm. Ranking and unranking techniques are useful for storage and retrieval of 
elements of S. Uniform random generation plays a role in Monte Carlo methods and in search 
algorithms such as hill climbing or genetic algorithms over classes of combinatorial objects. 
Uniform random generation of objects is always possible if one has an unranking algorithm 
since one can generate, uniformly, an integer in {0, • • • ,\S\ — 1} and unrank. 

We consider the set C n of trees with vertex set [n] = {1, . . . , n}. These trees are sometimes 
called Cayley trees and can be viewed as the set of spanning trees of the complete graph K n . 
Ranking and unranking algorithms for the set C n have been described by many authors. Indeed, 
efficient ranking and unranking algorithms have been given for classes of trees and forests that 
considerably generalize the Cayley trees (e.g., [3], [4], [5], [6], [7]). 

In a previous paper [14], we considered a more refined problem, namely, the problem of 
ranking and unranking subsets of C n with a specified degree sequences or a specified multiset of 
degrees. Let C n> i be the set of directed trees on V that are rooted at 1. That is, a directed tree 
T G C raj i has all its edges directed towards its root 1. We replace C n with the equivalent set 
U n> x. For any tree T £ C n , Yli=i ^ e 9T{i) = 2ra — 2. If s = (si, . . . , s n ) is a sequence of positive 



integers such that Y^h=1 Si = 2n — 2, then we let C n> g = {T£ C n) i : (degr(l), • • • 5 degxin)) = s}. 
Remmel and Williamson p3] proved that 

Similarly if S = {l ai , . . . , (n — l) Qn_1 } is a multiset such that Ya=i a « = 2ra — 2 and XIT=i a « = 
n, then we let C n> s = {T£ C nj i : {degxil), • • • , degrin)} = S}. It is easy to see from (llj) that 

. . . , a n y \si - 1, . . . , s n - 1J 

The basis of the ranking and unranking algorithms in [II] for C n g or C^, <j hinged on certian 
special properties of a bijection between C n> i and the class of functions J^n = {f ■ {2, . . . , n — 
1} —7- [n]} defined by Egecioglu and Remmel [3]. That is, in |14j . we proved that for any vertex 
i, 1 + equals that degree of i in the tree T = 6(/) when G(/) is considered as an 

undirected graph. This property allowed us to reduce the problem or ranking and unranking 
trees in C n ^g or C Hj s to the problem of ranking and unranking certain classes of set partitions 
of [n] . We were then able to modify known techniques for ranking and unranking set partitions 

HI] to construct efficient ranking and unranking algorithms for C n g or C n g. 

In this paper, we shall give efficient algorithms to rank and unrank two other natural subsets 
of C nt i, namely, the set of trees which have a given number of leaves or a prespecified set of 
leaves. If G = (V, E) is digraph and v £ V, we let 

indeg G (v) = \{u : (u,v) £ E}\, 

outdegc(v) = \{u : (v,u) £ E}\ and 

degc(v) = indegc{v) + outdegc{v). 

If T £ C n i, then we say that i is a leaf of T if and only if degxii) = 1- Fix k such that 
2 < k < n — 1 and C„ 1 equal the set of trees T in C n> i with k leaves. Similarly, if L is any 

subset of {1, ... , n} of size k, we let \ equal the set of trees T £ C n ,i such that i is a leaf of T 
if and only if i £ L. The main goal of this paper is to construct efficient ranking and unranking 
algorithms for the sets or C^ x . 

Just as in the case of the construction of the ranking and unranking algorithms for C n ^g 
or C n} s [13], the Egecioglu and Remmel bijection allows us to reduce the problem problem 
ranking and unranking algorithms or C\ 1 to the problem of finding ranking and unranking 
certain classes of set partitions. That is, we can restate the fundamental property of the bijection 
: T n -> C n ,i as 

deg QU) {i) = \f-\£)\ + l (3) 

for all / £ T n and i £ [n]. Now suppose that L is a subset of [n] of size k where 2 < k < n. Then 
by (E), it follows that if T £ and / = ®~ l (T), then f~ l {i) = if and only if * G L. Thus if 
J = [n] - L = {ji < . . . < j n - k }, then . . . , / _1 (jn-fc)) must be an ordered set partition 

of{2,...,n — 1} into n — k nonempty parts. Conversely, if we are given an ordered set partition 
7r = (it 1, . . . , 7r n _fc) of {2, . . . , n — 1} into n — k non-empty parts, we can define a function / £ T n 
such that = if and only if % £ L by setting f~ l {jt) = 7ft for t = 1, . . . , n — k. It follows 
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that \C, 'i \ is equal to the number of ordered set partitions of {2, . . . ,n — 1} into n — k parts. 
We then develop algorithms for ranking and unranking such classes of ordered set partitions by 
modifying ranking and unranking algorithms for the decreasing functions, permutations, and 
unordered set partitions found in [15j . 

Let S ni k denote the number of unordered set partitions of [n] into k parts. The numbers S nt k 
are call the Stirling numbers of the second kind and they satisfy the following recursion: 

S n ,k = if either k > n or n < 0, 
Sofi = 1, and 

Table 1 below gives the values of SVt,fc for 1 < n < 9. 
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Table 1 The values of SV^fc 
It follows from our arguments above that 



i a 



fc,Li 
n,l I 



(n - fc)!S' n _2 jn _jfc. 



(4) 



Similarly for any k such that 2 < k < n — 1, it is easy to see that C* 1 is the disjoint union of 



all Cn' L such that L C [n] of size k and hence 



(n - k)\S n - 2 ,n-k- 



(5) 



Note that both \C% \ and |C?*| can be as large as as 0{n\) so that the numbers involved in ranking 
and unranking can require 0{nlog{n)) bits. We show that our ranking and unranking algorithms 
require at most 0(n 2 ) comparisons of numbers y < n plus 0(n) operations of multiplication, 
division, addition, substraction and comparision on numbers x < \C n ' \ (x < |C^|). 
The outline of this paper is as follows. In Section 2, we describe the bijection : T n 



a 



n,l 

of [3] and discuss some of its key properties. In Section 3, we show that both and _1 can 
be computed in linear time. This result allows us to reduce the problem of efficiently ranking 
and unranking trees in C n ^g or C n ^s to the problem of efficiently ranking and unranking certain 
classes of ordered set partitions. In section 4, we shall recall the algorithms due to Williamson 
|15j for ranking and unranking decreasing functions, permutations and unordered set partitions 
which will be the building blocks of our final ranking and unranking algorithms. Finally in 



Section 5, we shall give our ranking and unranking algorithms for the sets C n 
examples. 



k,L 



or C„ and give 



3 



2 The G Bijection and its Properties 



In this section, we shall review the bijection : F n — > C n> i due to Egecioglu and Remmel (3] 



and give some of its properties. 

Let [n] = {1,2,..., n}. For each function / : {2, . . . , n — 1} — > [n], we associate a directed 
graph /, graph(f) = ([n],E) by setting E = {(i,f(i)) : % = 2, . . . , n — 1}. Following |13j . given 
any directed edge (i, j) where 1 < i, j < re, we define the weight of (i, j), W((i,j)), by 



where Pi,qi, Si,ti are variables for i = 1, . . . , n. We shall call a directed edge (i, j) a descent edge 
iii> j and an ascent edge if i < j. We then define the weight of any digraph G = ([n], E) by 



A moment's thought will convince one that, in general, the digraph corresponding to a function 
/ € F n will consists of 2 root-directed trees rooted at vertices 1 and n respectively, with all edges 
directed toward their roots, plus a number of directed cycles of length > 1. For each vertex v 
on a given cycle, there is possibly a root-directed tree attached to v with v as the root and all 
edges directed toward v. Note the fact that there are trees rooted at vertices 1 and n is due to 
the fact that these elements are not in the domain of /. Thus there can be no directed edges out 
of any of these vertices. We let the weight of /, W(f), be the weight of the digraph graph(f) 
associated with /. 

To define the bijection 0, we first imagine that the directed graph corresponding to / £ T 
is drawn so that 

(a) the trees rooted at n and 1 are drawn on the extreme left and the extreme right respectively 

with their edges directed upwards, 

(b) the cycles are drawn so that their vertices form a directed path on the line between re and 

1, with one back edge above the line, and the root-directed tree attached to any vertex on 
a cycle is drawn below the line between n and 1 with its edges directed upwards, 

(c) each cycle Cj is arranged so that its maximum element rrii is on the right, and 

(d) the cycles are arranged from left to right by decreasing maximal elements. 

Figure 1 pictures a function / drawn according to the rules (a)-(d) where n = 23. 

This given, suppose that the digraph of / is drawn as described above and the cycles of 
/ are ci(/), . . . ,c a (f), reading from left to right. We let ^ c ,(/) and l Ci (f) denote the right and 
left endpoints of the cycle Cj(/) for i = I,..., a. Note that if Cj(/) is a 1-cycle, then we 
let ?" Ci (/) = lci(f) be the element in the 1-cycle. ©(/) is obtained from / by simply deleting 
the back edges ( r c,(/)^Ci(/)) for i = l,...,a and adding the directed edges (?" Ci (/), l Ci+1 (f)) f° r 
i = 1, . . . , a— 1 plus the directed edges (n, l Cl (f)) and (r Ca ^, 1). That is, we remove all the back 
edges that are above the line, and then we connect n to the lefthand endpoint of the first cycle, 
the righthand endpoint of each cycle to the lefthand endpoint of the cycle following it, and we 
connect the righthand endpoint of the last cycle to 1. For example, ©(/) is pictured in Figure 2 




(6) 



w(G)= n w((i,j)). 



(7) 



4 



1 I gj ) 

2 15 

3 10 

4 20 

5 12 

6 11 

7 23 

8 22 

9 12 

10 3 

11 20 

12 14 

13 17 : 

14 5 , 

15 11 

16 3 • 

17 5 

18 1 

19 19 

20 6 

21 15 

22 3 




Figure 1: The digraph of a function. 
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Figure 2: 9(/). 



for the / given in Figure 1. If there are no cycles in /, then B(/) is simply the result of adding 
the directed edge (n, 1) to the digraph of /. 

To see that is a bijection, we shall describe how to define -1 . The key observation is 
that we need only recover that the directed edges {r c .if\, i Cj+1 m) for i = 1, . . . , a — 1. However 
it is easy to see that f ci m = vn\ is the largest element on the path from n to 1 in the tree ©(/)• 
That is, mi is then largest element in its cycle and by definition, it is larger than all the largest 
elements in any other cycle so that mi must be the largest interior element on the path from n 
to 1. Then by the same reasoning, r C2 ^ = mi is the largest element on the path from mi to 1, 
etc. Thus we can find mi, ... , m t . More formally, given a tree T G consider the path 

m = n, xi, . . .,mi,x 2 , . . . m 2 , . . . , x t , . . . , m t , 1 

where rrii is the maximum interior vertex on the path from m^-i to 1, 1 < i < t. If (mj_i, m^) is 
an edge on this path, then it is understood that Xj, . . . , mj = m^ consists of just one vertex and 
we define x« = wij. Note that by definition tuq = n > mi > . . . > m t . We obtain the digraph 
_1 (T) from T via the following procedure. 

Procedure for computing B X (T) : 

(1) First we declare that any edge e of T which is not an edge of the path from n to 1 is an edge 

of q-\t). 

(2) Next we remove all edges of the form (m t , 1) or (mj_i,Xj) for 1 < i < t. 
Finally for each i with 1 < i < t, we consider the subpath Xi, . . . , m,. 

(3) If mi = Xi, create a directed loop (mi, mi). 

(4) If mi / Xi, convert the subpath Xi, . . . ,m% into the 



5 



directed cycle Xi, . . . , rrii, X{. 



Next we consider two important properties of the bijection 0. First has an important 
weight preserving property. We claim that if ©(/) = T, then 

qnhW(f) = W(T). (8) 

That is, by our conventions, any backedge (fci(f)^ci(f)) are descent edges so that its weight is 
Qr itM f f \' Thus the total weight of the backedges is 

c i\J ) c i\J ) 

a 
i=l 

Our argument above shows that all the new edges that we add are also descent edges so that 
the weight of the new edges is 

a— 1 a 

QoKwtll ^(fA^hr^h = Qntll[qr c . (f) tl ai(f y (10) 
i=l i=l 

Since all the remaining edges have the same weight in both the digraph of / and in the digraph 
G(/), it follows that q n hW(f) = W{Q{f)) as claimed. 
It is easy to see that 

n-l 

W(f)=Y[[q i (t 1 + --- + t l )+ Pi ( Si+1 + --- + Sn )}. (11) 

f£T n i=2 

Thus we have the following result which is implicit in [3] and it explicit in |13j . 
Theorem 1. 

n-l 

W(T) = q n t 1 Y[[qi(ti + ---+ti)+pi(s i+1 + --- + s n )}. (12) 

Next we turn to a second key property of the O bijection. It is easy to see from Figures 
1 and 2 that deleting the back edges (?%(/)> 'c<(/)) for i = 1, . . . , o in graph(f) and adding the 
directed edges (^m, l Ci+1 (f)) for i = 1, . . . , a — 1 plus the directed edges (n, l Cl (f)) an d (fca(f)i 1) 
to get ©(/) does not change the indegree of any vertex except vertex 1. That is, 

indeg graph ( f) {i) = indeg® U) {i) for i = 2, . . . , n. (13) 

It is also easy to see that in going from graph(f) to ©(/), the indegree of vertex 1 increases by 
1, i.e., 

1 + indeg graph{f) (i) = indeg @(f) (l). (14) 

When we consider ©(/) as an undirected graph T, then it is easy to see that degxii) = 
outdegQtf\ + indegQ^y Thus since the outdegree of i in ©(/) is 1 if i ^ 1 and the outde- 
gree of 1 in ®(f) is zero, equations (13) and (14) imply the following theorem. 
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Theorem 2. Suppose that T is the undirected tree corresponding to G(/) where f € T n , then 
for i = 1, . . . , n, 

deg T (i) = l + \f- 1 (i)\. (15) 

ProofBy our definition of graph(f), it follows that indeg grap h(f){i) = \ f~ l {i)\ for i = 1, . . . , n. 
Thus by ( 13 ), for i = 2, . . . , n, 

deg T {i) = outdeg e{f) (i) + indeg e ^ f) (i) 
= 1 + indeg Q (fl(i) 
= l + indeg graph{f) (i) 

= i + ir l wi- 



Similarly by (14), 



deg T (l) = outdeg e{f) {l) + indeg B ( f) {l) 

= O + mde£/ (/)(l) 

= l + indeg gTaph{f) {l) 

= i + irHi)!- 



□ 



3 Construction of the spanning forests from the the function 
table in time 0(n) 

In this section, we shall briefly outline the proof that one can compute the bijections B and its 
inverse in linear time. Suppose we are given / £ T n . Our basic data structure for the function 
/ is a list of pairs (i, f(i)) for i = 2, . . . , n — 1. Our goal is to construct the directed graph of 
G(/) from our data structure for /, that is, for i = 1, ... ,n, we want to find the set of pairs, 
(i,U), such that there is directed edge from i to ti in 0(/). We shall prove the following. 

Theorem 3. We can compute the bijection O : T n — > C n> \ and its inverse in linear time. 

Proof. We shall not try to give the most efficient algorithm to construct ©(/) from /. In- 
stead, we shall give an outline the basic procedure which shows that one can construct ©(/) 
from / in linear time. For ease of presentation, we shall organize our procedure so that it makes 
four linear time passes through the basic data structure for / to produce the data structure for 

©(/)• 

Pass 1. Goal: Find, in linear time in n, a set of representatives ti,...,t r of the cycles of 
the directed graph of the function f. 

To help us find t\, . . . ,t r , we shall maintain an array A[2], A[3], ■ ■ ■ A[n — 1], where for each i, 
A[i] = (ci,pi, qi) is a triple of integers such Cj E {0, . . . , n — 1} and {pi, qi\ C {—1, 2, • • • , n — 1}. 
The Cj's will help us keep track of what loop we are in relative to the sequence of operations 
described below. Then our idea is to maintain, through the pi and qi, a doubly linked list of the 
locations i in A where q = 0, and we obtain pointers to the first and last elements of this doubly 
linked list. It is a standard exercise that these data structures can be maintained in linear time. 
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Initially, all the Cj's will be zero. In general, if Cj = 0, then pi will be the largest integer j 
such that 2 < j < i for which Cj = if there is such a j and pi = — 1 otherwise. Similarly, we 
set qi > i to be the smallest integer k such that n — 1 > k > i for which c& = if there is such 
a and qi = —1 if there is no such k. If a > 0, then q2 is the smallest integer j > 2 such that 
Cj = and qi = —1 if there is no such integer j. If c„_i > 0, then p n -\ is the largest integer 
k < n\ such that = and = — 1 if there is no such integer k. 

We initialize A by setting A[2] = (0, -1,92), ^[*] = (0, i - 1, i + 1) for m + 1 < i < n - 1, 
and A[n — 1] = (0,p n _i, —1). If 2 < n — 1 then (72 = 3 and p n -i = n — 2. Otherwise (2 = n — 1), 
these quantities are both —1. 

LOOP(l): Start with i x = 2, setting c 2 = 1. Compute /°(2), / 1 (2), / 2 (2), . . . , f kl (2), each 
time updating v4 by setting Cp( 2 ) = 1 an d adjusting pointers, until, prior to setting c^fcw 2 ) = 1) 
we discover that either 

(1) / fcl (2) g {l,n}, in which case we have reached a node in graph(f) which is not in the 
domain of / and we start over again with the 2 replaced by the smallest i for which a = 0, or 

(2) x = f kl (2) already satisfies c x = 1. This condition indicates that the value x has al- 
ready occurred in the sequence 2, /(2), / 2 (2), . . . , f kl (2). Then we set ti = f kl (2). 

LOOP (2): Start with %2 = q m +i which is the location of the first i such that Cj = 0, 
and repeat the calculation of LOOP1 with 12 instead of i\ = 2. In this manner, generate 
f°(i2), / 1 («2), / 2 (^2), • • • , / fc2 (*2), each time updating A by setting cp^ = 2 and adjusting 
pointers, until either 

(1) f kl {i2) £ {l,n}, in which case we have reached a node in graph(f) which is not in the 
domain of / and we start over again with the 12 replaced by the smallest i for which a = 0, or 

(2) x = f kl {i2) already satisfies c x = 2. (This condition indicates that the value x has al- 
ready occurred in the sequence 12, f (12) , f 2 («2 ) , • • • , f kl («2 ) • ) Then we set 1 2 = f kl {12 ) • 

We continue this process until q2 = — 1. At this point, we will have generated ti,...,t r , where the 
last loop was LOOP(r). The array A will be such that, for all 2 < % < n— 1, 1 < Cj < r identifies 
the LOOP in which that particular domain value i occurred in our computation described above. 

Pass 2. Goal: For i = 1, . . . ,r, find the largest element mi in the cycle determined by t{. 
It is easy to see that this computation can be done in linear time by one pass through the array 
A computed in Pass 1 above. At the end of Pass 2, we set /j = /(mj). Thus when we draw the 
cycle containing ti according to our definition of Oj(/), mi will be right most element in the 
and U will be the left most element of the cycle containing ti. However, at this point, we have 
not ordered the cycles appropriately. This ordering will be done in the next pass. 

Pass 3. Goal: Sort (Ii,mi),...,(I{.,mj) so that they are appropriately ordered according the 
criterion for the bijection ©(/) as described in by condition (a) -(d). 

Since we order the cycles from left to right according to decreasing maximal elements, it is 
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then easy to see that our desired ordering can be constructed via a lexicographic bucket sort. 
(See Williamson's book [15] for details on the fact that a lexicographic bucket sort can be carried 
out in linear time.) 

Pass 4. Goal: Construct the digraph of 0(/) from the digraph of f. 

We modify the table for / to produce the table for 0(/) as follows. Assume that (Zi, mi), . . . , (Z&, mk) 
is the sorted list coming out of Pass 3. Then we modify the table for / so that we add entries for 
the directed edges (n, Zi) and (m^, 1) and modify entries of the pairs starting with mi, . . . , 
so that their corresponding second elements are l 2 , ■ ■ ■ , Ik, 1 respectively. This can be done in 
linear time using our data structures. 

Next, consider the problem of computing the inverse of 0. Suppose that we are given the 
data structure of the tree T G C\, i.e. we are given a set of pairs , (i,U), such that there is a 
directed edge from i to t{ in T. Recall that the computation of _1 (T) consists of two basic steps. 

Step 1. Given a tree T £ C nj i, consider the path 

m = n,x 1 , . . . , m\,x 2 , ■ ■ ■ m 2 , . . . , x t , . . . , m t , 1 

where m« is the maximum interior vertex on the path from mj_i to 1, 1 < % < t. If (m^—i, m«) 
is an edge on this path, then it is understood that Xi, . . . ,rrii = rrii consists of just one vertex 
and we define Xi = mi. Note that by definition mo = n > m\ > ... > nit- 

First it is easy to see that by making one pass through the data structure for F, we can construct 
the directed path n —> a\ —)■... — >■ a r where 1 = a r . In fact, we can construct a doubly linked 
list (n, ai, . . . , a r _i, 1) with pointers to the first and last elements in linear time. If we traverse 
the list in reverse order, (1, a r _i, . . . , a±, n), then it easy to see that mt = a r _i, m tl is the next 
element in the list (a r -2) • • • , &i) which is greater than mt and, in general, having found m; = a s , 
then m.j_i is the first element in the list (a s _i, . . . , ai) which is greater than mj. Thus it is not 
difficult to see that we can use our doubly linked list to produce the factorization 

m = n, xx, . . . , m\,x 2 , ■ ■ ■ m 2 , . . . , x t , . . . , m t , 1 

in linear time. 

Step 2. We obtain the digraph -1 (T) from T via the following procedure. 
Procedure for computing 0J 1 (F) : 

(1) First we declare that any edge e of T which is not an edge of the path from n to 1 is an edge 
of -1 (T). 

(2) Next we remove all edges of the form (mt, 1) or (mj_i,Xj) for 1 < i < t. 
Finally for each % with 1 < i < t, we consider the subpath Xj, . . . , mj. 

(3) If mi = Xi, create a directed loop {mi, mi). 
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(4) If rrii Xi, then , convert the subpath Xj, . . . , into the 
directed cycle Xi, . . . ,rm, x%. 

Again it is easy to see that we can use the data structure for T, our doubly linked list, and 
our path factorization, tuq = n,x±, . . . ,mi,X2, ■ ■ -m?,, . . . ,xt, ■ ■ ■ ,mt,j to construct the data 
structure for graph(f) where / = _1 (T) in linear time. □. 

Given that we can carry out the bijection and its inverses in linear time, it follows that 
in linear time, we can reduce the problem of constructing ranking and unranking algorithms 
for C n to the problem of constructing ranking and unranking algorithms for the corresponding 
function class T n . 

4 Machinary for Ranking and Unranking Algorithms. 

Fix k such that 2 < k < n — 1. Recall that x equals the set of trees T in C rh \ with k leaves. 

Similarly, if L is any subset of {1, . . . , n} of size k, we let C*'i equal the set of trees T G C n> i 
such that i is a leaf of T if and only if i G L. The main goal of this section is to develop the 
basic machinary that is needed to give our final ranking and unranking algorithms for the sets 

Our ranking and unranking algorithms for x are based on six reductions. 

1. By the Egecioglu and Remmel bijection 6 of section 2, the problem of ranking and un- 
ranking trees in ± can be reduced to the problem of ranking and unranking the set F n ^ 
of functions / : 2, . . . , n — 1 — >■ [n] such that \{i G [n] : = 0}| = k. 

2. To specify the set of i G [re] such that f~ l {i) = for an / G F n j., we specify a decreasing 
function gf : {1, . . . , k} — > [n] whose range if {i : f~ l {i) = 0}. Let VF n) k denote the set of 
decreasing functions / : {1, . . . , k} — > {1, . . . , n}. Then given g G VT n ^, we let F„ iS equal 
the set of functions / : 2, . . . , n — 1 — > [n] such that {i G [n\ : = 0} = range(g). 

3. Given gf G T>T n ^, we specify a function / G F n>g , by giving an ordered set partition 7rj of 
[n— 2] into n — k parts. That is, suppose that tt = (iri, . . . , ir n -k) is an ordered set partition 
of [re — 2] into n — k parts and [n] — range(g) = {«!<■■■< i n -k}- Let ir + = (tt^, . . . ,^n-k) 
denote the set partition of {2, . . . , re — 1} which results from it by replacing each element 
x < n — 2byx + l. Then we can specify an / G F nj9 by declaring that f~ l {ij) = vr^ for 
j = 1,... ,n- k. 

4. To specify an ordered set partition of [n — 2] into n — k parts, we shall specify an unordered 
set partition V f = (Ti, . . . , T n _fe) of [re— 2] into n — k parts where miniVi) < • • • < miniVk) 
and a permutation Uj = a\ . . . (J n _k in the symmetric group S n _k- That is, given Tf and 
a j , we let T a be the ordered set partition of [re— 2] into k parts where T a = (T IJ1 , . . . , T (Tn _ k ). 

5. We shall associate to each permutation a = a\ . . . a n G S n , a sequence, h a , called the direct 
insertion sequence associated with a. We defined h a = • • • , h(n)) by recursion as 
follows. First if a G Si, h a = (0). If re > and a = a\ ...,£>"„ G S n where o\ = j, then let 
o~ be the permutation derived from o<i . . .a n where we replace each i > j in the sequence 
(T2 ... <r n by i — 1. For example, ifcr = 41256 3, then o~ = 1 2 4 5 3. This given, 
suppose /i CT - = (/i~(l), . . . ,h~(n — 1)), then we define h a = (j — 1, h~(l), . . . , h~(n — 1)). 
For example, <r = 4 1 2 5 6 3, then h a = (3, 0, 0, 1, 1, 0). 
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6. We shall associate to each unordered set partition T = (Ti, . . . , r n _&) of [n — 2} into n — k 
parts where min(Ti) < • • • < min(T n _k), a sequence sr = • ■ • , s( n — 2)) associated 

to its corresponding restricted growth function. That is, s(i) = j — 1 if and only if i G Tj. 
For example, suppose T = ({1, 4}, {2, 3, 8}, {5}, {6, 7}), then s 7 = (0,1,1,0,2,3,3,1). 

It follows that we can specify any / G F n> k by a triple (g, a, T) where g GG VT^k, cr G S n ~k 
and T is an unordered set partition of [n — 2] into n — k parts. We can thus identity / with a 
triple of sequences 

Seq(f) = <( 5 (1), . . . , g(k)), (h(l), ...,h(n-k)), ( a (l), ...,s(n- 2))) 

where /i CT = (h(l), . . . , /i(n — fc)) and sr = (s(l), • • • , s(n — 2)). We can then order the functions 
/ £ Fn,k according to the lexicographic order of their associated sequences Seq(f). 

9 4 7 

T = . a* . 

• 

S 




1 (T) = 




8 3 5 6 



Figure 3: A tree in Cg X . 

For example, consider the tree in T G Cg 1 pictured in Figure [3J We have also pictured 
_1 (T) = / so that / can be specified by its table of preimages / _1 (1), . . . , /~ 1 (9) given in 
Table 2. 



r x 0) 


r i (2) 


r x (3) 


r A (4) 


r A (5) 


/ _1 (6) 


/ _1 (7) 




r x (9) 


{5,6} 


{2,3} 





{7} 








{4,8} 









Table 2 



Thus / is associated with the decreasing functions gj G VTg^ given by gf = (9,8,6,5,3) 
and the ordered set partition irj = ({5, 6}, {2, 3}, {7}, {4, 8}) of {2,..., 8} into 4 parts. Let 
Hf = ({4, 5}, {1, 2}, {6}, {3, 7}) be the ordered set partition of [7] into 4 parts which results 
by replacing each element i in tt^ by i — 1. Then 7rj is specified by the underlying unordered 
set partition Tf = ({1, 2}, {3, 7}, {4, 5}, {6}) and the permutation aj = 3 1 4 2. Finally cry 
associated to insertion order sequence h a . = (2,0,1,0) and T is associated to the restricted 
growth function sp, = (0, 0, 1, 2, 2, 3, 1). Thus 

5eg(/) = ((9, 8, 6, 5, 3), (2, 0, 1, 0), (0, 0, 1, 2, 2, 3, 1)). 

In the next subsection, we shall provide the basic lemmas about ranking and unranking leaves 
of trees which will allow us to reduce the problem of ranking Seq(F n> k) = {Seq(f) : / G -F n ,fc} 
according to lexicographic order to the problems of ranking and unranking decreasing functions 
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according to lexicographic order, of ranking and unranking insertion sequences of permutations 
according to lexicographic order, and of ranking and unranking restricted growth functions 
according to lexicographic order. 



4.1 Basic Lemmas for Ranking and Unranking Algorithms. 

To develop our ranking and unranking algorithms for C„ 1 and C n \ , we first need to make some 
general remarks about ranking and unranking paths in planar trees. Given a rooted planar tree 
T, let L(T) be the numbers of leaves of T and Path(T) be the set of paths which go from the 
root to a leaf. Then for any path p £ Path(T), we define the rank of p relative to T, rankrip), 
to be the number of leaves of T that lie to the left of p. 

Given two rooted planar trees T\ and T2 , Ti (g T2 is the tree that results from T\ by replacing 
each leaf of T\ by a copy of T2, see Figure 3. If the vertices of T\ and T2 are labeled, then 
we shall label the vertices of T\ (g T2 according to the convention that each vertex v in T\ 
have the same label in T\ <g T2 that it has in T\ and each vertex w in a copy of T2 that is 
decendent from a leaf labeled / in T\ has a label (l,s) where s is the label of w in T2. Given 
rooted planar trees T±, . . . , T& where k > 3, we can define T± (g T2 (g ■ ■ ■ (g by induction as 
(Ti (g) ■ ■ ■ (g) Tfc_i) (g) Tfc. Similarly if T\ . . . , are labeled rooted planar trees, we can define the 
labeling of T\ ® T2 ® • • • ® by the same inductive process. 




Figure 4: The operation Xi <g> T 2 . 



Now suppose that we are given two rooted planar trees T\ and T2 and suppose that p\ G 
Path(T\) and P2 £ Path(T2). Then we define the path pi <g) ^2 in Ti <g> T2 which follows p\ to 
its leaf I in Ti and then follows P2 in the copy of T2 that sits below leaf / to a leaf (I, I') in 
T\ (g T2. Similarly, given paths 6 Tj for i = 1, . . . k, we can define a path p = p± (g ■ ■ ■ <g p^ € 
Path{T\ (g T2 (g ■ ■ ■ (g Tfc) by induction as (pi (g • • • (g Pfc-i) <g Pfe- 

Next we give two simple lemmas that tell us how to rank and unrank the set of paths in 
such trees. 

Lemma 4. Suppose that T±, . . . , T& are rooted planar trees and T = T\ (g T2 (g ■ ■ ■ (g T^. T/ien 
/or any pai/i p = pi (g • • • (g Pfc G Path{T), 

k k 

rank T (p) = Y,rank Tl {p j ) ]J L(T,) (16) 
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Proof. We proceed by induction on k. Let us assume that T±, . . .Tj. are labeled rooted planar 
trees. 

First suppose that k = 2 and that p\ is a path that goes from the root of T\ to a leaf labeled 
li and p2 goes from the root of P2 to a leaf labeled li. Thus p\ <g> P2 goes from the root of 
T\ (g) T2 to the leaf l\ in Ti and then proceeds to the leaf (Ji, l 2 ) in T± (g) I2. Now for each leaf 
/' to the left of l\ in T2, there are L{P2) leaves of T\ <g> P2 that lie to left of (Ji, I2) coming from 
the leaves of the copy of T2 that sits below I'. Thus there are a total of L(T2) ■ rank^ipi) such 
leaves. The only other leaves of P\ T2 that lie to left of p\ <%> P2 are the leaves of the form 
(h,l") where I" is to left of P2 in T2. There are rankT 2 {P2) such leaves. Thus there are a total 
of rankT 2 {P2) + ^(^2) ■ rank^ipi) leaves to left of pi ®p2 and hence 



rankr^TiiPi ® P2) = rank T2 {p2) + L(T 2 ) ■ rank^ipi) 



as desired. 

Next assume that (16) holds for k <n and that n > 3. Then 

rank Tl ®-®T n {Pi ®---®Pn) = 
rank {Tl ®...® Tn _ l) ® Tn (p l <g) • • • <8>p n -\) ® p n ) = 

n— 1 n— 1 

rank Tn (p n ) + L(T n )C^2rank T] (pj) ] [ L(T/)) = 

i=l J=j+i 

n n 

Y^rank Tj {pj) JJ L(Tj). 
i=i Z=j+1 

□. 

This given, it is easy to develop an algorithm for unranking in a product of trees. The proof 
of this lemma can be found in 



Lemma 5. Suppose that T\, . . . , Tj- are rooted planar trees and T = T\ (g T2 <g • • • <g T k . Then 
given ap £ PathiT) such that rankx(p) = ro, p = j?i<g- ■ -®Pk £ Path(T) where rankx^Pi) = q% 
and 



qxY[L(Ti) + n where < r x < IIL2 L ( T 0> (17) 

Z=2 
fc 

ri = 92 [J + r 2 w/jere < r 2 < fits L (Ti), (18) 

Z=3 

; (19) 

Tk-2 = qk-\L{T k ) + r fe _i w/iere < r fc _i < L(T fc ) and (20) 
rk-i = qk (21) 

Our idea is to construct trees Tx>T n k 1 Txo n and Tjig n k so that 

1. the paths of TvT n k correspond to the decreasing functions in P>T n ^ ranked according to 
the lexicographic order, 
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2. the paths of Txo n correspond to the permutations of S n ranked according to the lexico- 
graphic order on their insertion sequences, and 

3. the paths of Tjig n k correspond to set restricted growth functions TZG n ,k ranked according 
to the lexicographic order. 

It will then easily follow that the paths of 

T<DT ntk ® T XOk ® T n g n _ 2 k 

naturally correspond to the sequences in Seq{F n ^) = {Seq(f) : / G F nk } ranked according 
to lexicographic order. Thus we can use Lemmas [4] and [5] to obtain a ranking and unranking 
algorithms to Seq{F nk ) relative to the lexicographic order once we have constructed the trees 
T-DT n k ; TxOn an d Tng n k and developed ranking and unranking algorithms for them. Our 
next three subsections will be devoted to constructing the trees Tt>T n k , 3xo n and Tjig n k and 
specifying ranking and unranking algorithms for them. 

4.2 Ranking and Unranking Decreasing Functions. 

In this subsection, we consider the problem of ranking and unranking the set T>J- nt k of decreasing 
functions / : {l,...,k} — > {1, . . . ,n} relative to lexicographic order. A number of authors 
have developed ranking and unranking algorithms for T>J- n ^- We shall follow the method of 
Williamson |15| . First, we identify a function / : {1, ... , k} — > {1, ... ,n} with the decreasing 
sequence (/(l), . . • , f(k)) where n > /(l) > . . . > f(k) > 1. We can then think of the sequences 
as specifying a path in a planar tree Tx>T n k which can be constructed recursively as follows. At 
level 1, the nodes of Tx>j? n k are labeled k, . . . ,n from left to right specifying the choices for /(l). 
Next below a node j at level one, we attach a tree corresponding to Tx>p j _ 1 k _ 1 where a tree 
TvTa i consists of a tree with a single vertex labeled a. Figure 3 pictures the tree Tt>f 63 . 




Figure 5: The tree Tt>f 73 - 

Then the decreasing sequence (6,2,1) corresponds to the path from the root to the node 
which is specified with an arrow. It is clear that the sequences corresponding to the paths the 
tree Tx>j^ 63 appear in lexicographic order from left to right. Thus the rank of any sequence 
(/(l), . . . , f(k)} E VFn^u is the number of leaves of the tree to the left of the path corresponding 
to (/(l), • • • , f(k)). Hence the sequence (6, 2, 1) has rank 10 in the tree Tpjr e 3 . 

This given, suppose we are given a sequence (/(l), • • ■ , /(&)) in Tpjr n k . Then the number of 

leaves in the subtrees corresponding the nodes k, . . . , /(l)— 1 are respectively (jjjljO , L^i) > • • • , 
Thus the total number of leaves in those subtrees is 

(l:m*-0+-+( / ivm / t i ) 
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Here we have used the well known identity that Yus=k-i (fc-i) = (D - ^ follows that the rank of 
(/(l), . . . , f(k)) in T VTn k equals plus the rank of (/(2), . . . , f(k)} in T^ /(1) _ l fc _ r The 

following result, stated in [To], then easily follows by induction. 

Theorem 6. Let f : {1, . . . , k} — > {1, . . . , re} &e a descreasing function. Then the rank of f 
relative to the lexicographic order on T)T n & is 

™„^(^(«l>(^- 1 ) + ... + (^- 1 )- 

It is then easy to see from Theorem [6j that the following procedure, as described by 
Williamson in [15], gives the unranking procedure for VJ- n ^. 

Theorem 7. The following procedure U N RAN K (m) computes 

f = (f(l),...,f(k)) 
such that Rank-£)F n k (f) = m for any 1 < k < n and < m < (£) - 1. 

Procedure UNRANK(m) 

initialize m' := m, t := 1, s := k; (1 < k < n, < m < (^) — 1) 
while £ < A; do 

begin f(t) — 1 = max{y : < m'}; 

m':=m'- 
i := t + 1; 
s := s — 1; 

end 

4.3 Ranking and Unranking Permutations 

The problem of ranking and unranking permutations according to lexicographic order on inser- 
tion sequences is quite easy. We can then think of the insertion sequence h a = (/i(l), . . . , h(n)) 
of a permutation a G S n as specifying a path in a planar tree Tjc^ which can be constructed 
recursively as follows. At level 1, the nodes of Txo n ar e labeled 0, ... ,n — 1 from left to right 
specifying the choices for h\. Next below a node j at level one, we attach a tree corresponding 
to Txo n -i where a tree Tjo 1 consists of a tree with a single vertex labeled 0. Figure [6] pictures 
the tree Ts 3 - 

This given, it easy to see from our definitions that the following theorems, which can be 
found in [15J, hold. 

Theorem 8. Let S n denote the symmetric group of all permutations of [n] . Order S n by defining 
a < t iff h a <i ex h T , then the rank of a, Rank(a) = \{t £ S n : r < a}\, can be computed by the 
formula 

n-l 

Rank[a) = ^h(i) ■ (n-i)\ (23) 

k=l 

ifh a = (h(l),...,h(n)). 
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o 4 o4 i)4 o4 i)4 o4 

Figure 6: The tree Tjo 3 - 

Theorem 9. Let S n denote the symmetric group of all permutations of [n] . Order S n by defining 
a < t iff h a <iex h T . Then if < p < n\ — 1, then we construct the insertion sequence 
h a = (h(l), . . . , h(n)) of the permuation a G S n such that Rank(a) = p as follows. Let p = po, 
then set h n = and compute the following k = 0, . . . , n — 2: 

Po = h(l)(n — 1)! + p\ where p\ < (n — 1)! 
pi = h(2)(n — 2)1 + p2 where p2 < (n — 2)1 

Pk = h(k + l){n — k)\ + pk where pk < (n — k)\ 



5 Ranking and Unranking Algorithms for Unordered Set Par- 
titions. 

The main purpose of this section is to give ranking and unranking algorithms for the set of 
all unordered set partitions of [n] into k parts, S Uj k- We shall follow the general approach of 
Williamson [15] and give algorithms to rank and unrank the set of surjective restricted growth 
functions relative to lexicographic order. A restricted growth function / : [n] — > {0, 1, . . . , k — 1} 
is a function such that 

(i) /(l) = and 

(ii) for all 1< i < n, f(i) < 1 + moz({/(l), ...,/(*- 1)})- 

We let TZQn^k denote the set of all restricted growth functions / that map [n] onto {0, 1, . . . , k — 
1}. There is a natural one-to-one correspondence I between TZQ n ^ and S n ^. That is, suppose 
7r = {tti, . . . ,TTk) is a set partition of [n] into k parts where min(iri) < . . . < min^k)- Then 
define L(jr) = f w 6 TZG n ,k by letting f n (i) = j — 1 where i G vr^ for all i G [n\. It is easy to see 
that if / G Ug njk , then = tt/ = (/ _1 (0)> • • • > f~ l { k ~ We shall identify a restricted 

growth function / G TZQ n ,k with the sequence (/(l), • . . ,f(n)} and then order lZQ n ^ via the 
lexicographic order on such sequences. 

Let (/(l), . . . , f(n — j)), < j < n — 1, be a restricted growth function and suppose that 
max({f(l), . . . , f(n — j)}) = m where < m < k — 1. Let 

£«(j,m,/(l),...,/(n-j))= (24) 
{(x n _j + i,. . . ,x n ) | (/(!),..., /(n - j),x„_ i+ i, . . . ,x n ) G UQ n> k}- 
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The set £ ( k \j, m, /(l), ■ • • , f(n — j)) represents all ways of extending the sequence 
/(l), . . . , f(n — j) to a surjective restricted growth function. Note that 

(1) f( fc )(0,m,/(l),...,/(n)) = 0for m^k-1, 

(2) £( k \0, k — 1, /(l), • • • , /(«)) = {e} where e is the empty string, and 

(3) for < j < n - 1, 

£«(j,m,/(l),...,/(n-j)) = 

m 

\Jt*£^(j-l,m,f(l),...,f(n-j),t) 
t=o 

U(m + 1) * - 1, m + 1, /(l), . . . , f(n - j),m + 1). 

This union is disjoint. The notation t * X means concatenate t with each element in X. Recall 
that we have assumed above that mox({/(l), . . . , f(n — j)}) = m. In the sets 

£^(j-l,mj(l),...,f(n-j),t) 

we have < t < m, thus in the function 

</(l), ...J(n-j), f(n - j + 1)) = </(l), ...,f{n- j),t) 

we still have max({f(l), . . . , f(n — j + 1)}) = m. Given that max({f(l), . . . , f(n — j)}) = m, 
we can also assign f(n — j + l) = m + l and still have (/(l), . . . , f{n — j + 1)) a restricted growth 
function. The set (m + 1) * £( k \j — 1, m + 1, /(l), . . . , f(n — j),m+ 1) includes all such elements 
in £( k \j, m, /(l), . . . , /(n — j)). Note that £^(j, m, /(l), . . . , /(n — j)) does not depend on the 
actual values of /(l), . . . , f(n — j). Setting 

\£^ k \j,mJ(l),...,f(n-m=E {k Hj,m) 

we have the following recursion 

(j, m ) = (m + (j - 1, m) + ,E (fc) (j - 1, m + 1) 

The initial values are E(0, m) = if m 7^ k — 1, i£(0, fc — 1) = 1. This recursion will be the basis 
for our ranking and unranking procedures for TZG n ,k and hence for 5„ 5 fc. 
For example, Table 3 below lists the values of E^\m,n). 



m 

n 





1 


2 


3 


4 


5 


6 


7 














1 














1 








1 


4 













2 





1 


7 


16 












3 


1 


9 


37 


64 











4 


10 


55 


175 


256 










5 


65 


285 


781 












6 


350 


1351 














7 


1701 

















Table 3 E^ A \m,n) 
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It is easy to construct a planar tree T n whose paths correspond to restricted growth functions 
(a n -i,a n -2, ■ ■ ■ , ao) by having the root of the T n labelled with and by induction if a node 
corresponds to a path (a n _i, . . . , at+i), then the descendants of that node would be labeled from 
left to right with 0, . . . , mt, mt+1 where mt = max({a n -i, . . . , at+i}); see Figure reffigure:RecT. 
The tree T-jig n k would consists of those nodes which lies on a path (a ra _i, a n -2> ■ • ■ > «o) such 
that max({a n -i, a n -2, • • • , ao}) = k—1. For example, = 6 and the tree Tng 4 3 is pictured 
in Figure [7] 




Figure 7: The tree T n g i3 . 

More generally, it is easy to see that if we are at node corresponding to the sequence 
(a n _i, . . . , at+i) in Tjig n k where m t = max(a n -i, . . . , ott+i) < k — 1, then, as pictured in figure 
[8j our choices for at is either 0,1,... ,mt,mt + 1. Moreover by our definition of E^(t,mt), 
the number of leaves in each of the subtrees where we branch to either 0, 1, . . . ,mt is exactly 
(t, mt) and the number of leaves in the subtree where we branch to mt + 1 is (t, mt + 1). 
It follows that if we take the branch corresponding to j out the node corresponding to the se- 
quence (a n _i, . . . ,at+i), then the number of nodes in the subtrees to the left of that branch 
is j ■ E^ k \t, mt). If we are at node corresponding to the sequence (a n _i, . . . , at+i) in T-jig n k 
where mt = max(a n -i, . . . , at+i) = k — 1, then we can only branch to either 0, 1, . . . , mt but 
it is still the case that, if we take the branch corresponding to j out the node correspond- 
ing to the sequence (a n _i, . . . ,at+i), then the number of nodes in the subtrees to the left of 
that branch is j ■ E^ k \t,mt). The following result of Williamson [JS] for ranking sequences 
/ = (a„_i,a n _2,...,ao) be an element of TZG n ,k relative to lexicographic order, then easily 
follows by induction. 




E(t,m t ) E(t,m t ) E(t,m t ) E(t, m t + 1) 



Figure 8: Typical node in T-jig n k . 
Theorem 10. Let f = (a n -i, ««-2, • • ■ , «o) be an element oflZQ n ^ and for each t = 0, 1, ... , n— 
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2, let nit = max({a n -i, . . . ,at+i})- Then 

n-2 

rank n g n k (f) = £ a t £« (i, m t ), (25) 
t=o 

For example, suppose that / = (0,1,0,2,0,3,1,2) = (aj,...,ao) G TZGs^- Then for i = 

0, ...,6, we let wi; = max({ci!7, . . . , c^+i}) so that mo = 3, mi = 3, 777,2 = 2,7773 = 2,m4 = 

1, m5 = 1 an d m6 = 0. Thus using the values of E^(n,m) from Table 3, we see that 

6 

rank nSnk (f) = ^a 4 ^ (4) (t,m t ) 
t=o 

= 2 • E^(0, 3) + 1 • E<-*>(1, 3) + 3 • £ (4) (2, 2) + • £ (4) (3, 2) + 

2 • E^(4, 1) + • £ (4) (5, 1) + 1 • £ (4) (6,0) 
= 2-1 + 1-4 + 3-7 + 0-37 + 

2 • 55 + • 285 + 1 • 350 
= 487. 

One can also use Figure [8] to give an inductive proof of the validity of the following unranking 
algorithm for the element of TZG n ^ relative to lexicographic order. 

Theorem 11. Let < r < S n k where S n k is the number of set partitions of [n] into k parts. 
Then we can construct f = (a n _i, a n -2, ■ ■ ■ , ao) G TtGn.k such that rankng n k (/) = r as follows. 

Step 77 — 1 Set a n -i = 0. 

Step n-2 Set m n _2 = and let a n _2 equal the largest s such that < s < m n _i + 1 and 

< r. Then set 

r : r - a n _ 2 -E^(n ~ 2 ,"in-i)- 

Step t Assume that we have defined a n -i, . . . , at+i- Let mt = max({a n -\, . . . ,ctt+i}). Then 
set at equal the largest s such that < s < mt + 1 and sE^ k '(t,mt) < r. Then set 

r : r — a t E^ k \t, m t ). 

Thus, for example, to compute the / = (a 7 , . . . , ao) G TZGsa whose rank is 1000, we would 
carry out the following steps. 

Step 7. We set r = 1000 and a 7 = 0. 

Step 6. m 6 = a 7 = 0. E^(6,0) = 350. Since 1 • E^(6,0) = 350 < 1000, then a 6 = 1 and we 
set r = 1000 - 350 = 650. 

Step 5. m 5 = max({0, 1}) = 1. E^(5, 1) = 285. Since 2 • E^(5, 1) = 570 < 650, then a 5 = 2 
and we set r = 650 — 570 = 80. 

Step 4. 7774 = max ({0,1,2}) = 2. £( 4 )(4,2) = 175. Since • £( 4 '(4,2) = < 80 < 1 • 
£( 4 )(4, 2) = 175, then a 4 = and we set r = 80 - • 175 = 80. 
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Step 3. m 3 = max({0, 1,2,0}) = 2. #( 4 )(3,2) = 37. Since 2 • #( 4 )(3,2) = 74 < 80 < 3 • 
£W(3, 2) = 111, then a 4 = 2 and we set r = 80 - 2 • 37 = 6. 

Step 2. m 2 = max({0,l,2,0,2}) = 2. £< 4 )(2,2) = 7. Since • E^(2,2) = < 6 < 1 • 
E^(3, 2) = 7, then a 2 = and we set r = 6 - • 7 = 6. 

Step 1. mi = max({0, 1,2,0,2,0}) = 2. £( 4 )(1,2) = 1. Since 3-^ 4 )(l,2) = 3 < 6, then 
a\ = 3 and we set r = 6 — 3 • 1 = 3. 

Step 0. m = max({0, 1,2,0,2,0,3}) = 3. E^(0,3) = 1. Since 3 ■ E^(0,3) = 3 < 3 < 
4 • E^(0, 3) = 4, then a = 3 and we set r = 3 - 3 • 1 = 0. 

Thus the restricted growth function / = (0,1,2, 0, 2, 0, 3, 3) is the element of rank 1000 in TZGs,4- 

6 Ranking and Unranking Algorithms for C^'f or 

The main goal of this section is to give our final ranking and unranking algorithms for the sets 
We start by considering the ranking and unranking algorithms for 1 . Recall 

\^n,l\ = Q X (n-k)\ X ,S„_2,n-fe- 

Ranking Algorithm for C*i- 

Step 1 Given T G C^, construct / = O'^T). 

We assume that we start with the data structure for a tree T 6 x which consists of pairs 
(i, j) such that (i, j) is a directed edge in T. We then compute / = 6 _1 (T). By our results in 
section 3, we can construct the set of pairs such that f(i) = j in linear time. 

Step 2. Find g G VT n , k , h G XO„_ fc , and s G 1ZQ n - 2 ,n-k such that Seg(/) = ((y(l), . . . g(k)), (h(l), 
fc)),(s(l),..., S (n-2))). 

Note that by making one pass through the data structure for /, we can construct of table 
Table(f) = ((1, Vi), . . . , (n, V n )) where for each i, V t = if = or ^ a linked list of 

the elements of / _1 (i) in increasing order if / _1 (i) / 0. For example if we start with the tree 
T G (7g ;1 and / = _1 (T) pictured in Figure 3, then one can see from Table 1 that we would 
produce the following table. 

Table(f) = ((1, (5, 6)), (2, (2, 3)), (3, 0), (4, (7)), (5, 0), (6, 0), (7, (4, 8)), (8, 0), (9, 0)). (26) 

We can then read Table(f) from right to left to construct the decreasing sequence 

(g(l),...,g(k)) 

of those i such that f^ 1 (i) = 0. Such elements correspond to the leaves of the tree T by Theorem 
2. Thus is clear that we can construct (g(l), . ■ ■ ,g(k)) from T in linear time. For the tree T 
pictured in Figure 3, we would produce 

( 5 (1),..., 5 (5)) = (9,8,6,5,3). (27) 
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By reading Table(f) from left to right, we can construct the ordered set partition of {2, . . . , n— 
1} into n — k parts, ir + (f) = (vrj 1 ", . . . , vr^"_ fc ), which consists of the V^'s which are nonempty in 
Table(f) and a set of pointer from 7r+ to j. For our example, 

tt + (/) = ((5,6),(2,3),(7),(4,8)). (28) 
Next we make a pass through ir + (f) and create 

(i) a new ordered set partition ir(f) = (tti, . . . ,ir n -k) of [n — 2] into n — k parts along with 

pointers from TTj to j by replacing each element i by i — 1, 

(ii) a sequence u(f) = (u(l), . . . , u(n — k) where u(i) = min(irf) — 1, and 

(iii) a sequence v(f) = (v(l), . . . ,v(n — 2)) where v(i) = if i appears in u(f) and v(i) = 1 
otherwise. 

For our given example of /, we would produce 

tt(/) = ((4,5),(1,2),(6),(3,7)), (29) 
«(/) = (4,1,6,3), (30) 

and 

v(f) = (0,1,0,0,1,0,1). (31) 

Again, it is easy to see that we can construct 7r(/), u(f), and v{f) in linear time from T. 

Next we use v (/) to construct 5(f) = (5(1), . . . , 5(n— k)) where 5(1) = and for 1 < i < n— 2, 
<K*) = S}=i Thus 5(i) is the number of j < i such j is not a minimum element in one of 

the parts ir(f). It follows that if we let a = (cr(l), . . . ,a(n — k)) = (u(l) — 5(u(i)), . . . ,u(n — 
k) — 5(n — k)), then a is a permutation of n — k which represent the relative order of parts of 
7r(/) according to increasing minimal elements. For example, for our given example of /, 

<5(/) = (0,0,l,l,l,2,2) (32) 

and 

a = (4 - 1, 1 - 0, 6 - 2, 3 - 1) = (3, 1, 4, 2) (33) 

Note that if we order the parts of vr(/) according to increasing minimal elements to get n~(f) = 
(7rf , . . . , ir~_ k ), the 7Tj = 7r~^y For our given example of /, 

7r-(/) = ((l,2),(4,5),(6),(3,7)). (34) 

Note that we change the pointer on 7r(/) = (iri, . . . , TT n -k) so that TTj pointers to Oj — 1, the part 
of 7r(/) which has the smallest minimal element will point to 0, the part of 7r(/) which has the 
second smallest minimal element will point to 1, etc. It follows that if we make another pass 
through 7r(/) and as we encounter an element x, we record Oj — 1 where the part that contains 
x in 7r(/) points to j, we can construct the sequence s(f) = (s(l), . . . , s(n — k)), then s(f) will 
be the restricted growth function associated to the unorderd set partition ir~(f). For our given 
/, this process would produce 

s(f) = (0,0,1,2,2,3,1). (35) 
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Again it is easy to see that we can produce s(f) in linear time from T. 

Finally we need to construct the direct insertion sequence h(f) = (h(l), . . . , h(n — k)) for a. 
We can produce the insertion sequence by recursion. That is, h(l) = <r(l) — 1 and (h(2), . . . , h(n— 
k)) is the insertion sequence for the permutation r = (r(l), . . . , r(n — k — 1)) where r(i) = a(i+l) 
if a(i + 1) < cr(l) and r(i) = a(i + 1) — 1 if a(i + 1) > c(l). Thus it requires 0((n - k) 2 ) steps 
to produce h(f). For our given example of /, 

h(f) = (2,0,1,0). (36) 

It follows that it require 0(n 2 ) step to produce Seq(f). 



Step 3. Use the algorithms of section 4 to compute rankx>T n k (g(f)), rankzo n _ k (h(f)), and 
rank n g n _ 2n _ k {s{f)). 

It is easy to see from our ranking algorithms of section 4 that if we start with table of binomial 
coefficients and E^ n ~ k \r, s), then it requires 0(n) operations of addition, multiplication, and 
comparisions to compute each of rank VJ r n k (g(f)), rank X o„_ k (h(f)), and rankng n _ 2tn _ k (s(f)). 
Since the integers involved have length at most 0(nlog(n)), it follows that it requires 0(n 2 log(n)) 
steps to compute rank VJ r nk (g(f)), rank XOn _ k (h(f)), and rank n g n _ 2n _ k (s(f)). 
For our given /, 

rank v ^Mf)) = £ (^jA ( 37 ) 



i=i 



= 56 + 35 + 10 + 6 + 2 = 109, 

4 



rankxoMf)) = £>«(«-*)! (38) 

i=l 

= (2 x 3!) + (0 x 2!) + (1 x 1!) + (0 x 0!) 
= 7, 

and if s(f) = (a^, a$, . . . , oto) and m,i = max({a.Q, . . . , Oj+i}), then 

5 

rank n g 7A = £ a^ 4 ) (i, m,) (39) 

i=0 

= (0 x £ (4) (5,0)) + (1 x £ (4) (4,0)) + (2 x £ (4) (3,1)) 

+(2 x E&{2,2)) + (3 x £ (4) (1,2)) + (1 x ^(0,3)) 

= (0 x 65) + (1 x 10) + (2 x 9) + (2 x 7) + (3 x 1) + (1 x 1) = 46. 



Step 4 It then follows from our analysis in section 4 that 

rank & = rank V T n k {g{f)) ■ (n - k)\S n - 2 , n -k (40) 

71,1 

+rank n g n _ 2n _ k {s(f)). 
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For our given /, n = 9, k = 5, (n - k)\ = 4! = 24, S 7 , 4 = 350, 4! x 350 = 8400. Thus 

rank^ (T) = (109 x 8400) + (7 x 350) + (46) = 918, 096 (41) 

°9,1 

Thus the T pictured in Figure 3 is the tree with rank 918,096 in Cf^ 

The unranking algorithm for 1 is relatively straight forward given the unranking algo- 
rithms discussed in section 4. 



Unranking Algorithm for C\ x . 

Problem: Given r such that < r < (£) x (n — k)\ x S n -2,n-k, find T G C"^ such that 
rank^ k (T) = r. 

Step I Find r\ r 2 , and r3 such that 

r = n x (n - /c)!S'„_2,„-fc + ui where < ni < (n - fc)!S , „_ 2 ,n-fe 

= r 2 x S' n _2,n-fc + ^2 where < u 2 < S n _ 2 ,n-k 
r 3 = u 2 - (42) 

Step II Find g G VT n ,k, h G XO n -k and s G HQ n -i,n-k such that 

n = rank VTnk {g) (43) 

r 2 = ranfc I0n _ t (/j) (44) 

r 3 = ran^„_ 2n _ fc (s) (45) 

Step III Find / G F n such that 

Seq(f) = {(g(l), g(k)), (h(l), ...,h(n-k)), (s(l), ...,s(n- 2))). (46) 

Step IV Let T = 9(/). 

Note that since we are dealing with numbers whose length is of 0(nlog(n), Step I requires 
0(nlog(n)) steps. The unranking algorithms for DJ^n^, TO n _k and TtG n ~2,n-k each require 
0(n) operations of addition, mulitiplication, division and comparisons. Again, since we are 
dealing with numbers whose length is of 0(nlog(n), Step II requires 0(n 2 log(n)) steps. It is not 
difficult to see that we can reconstruction / from g, h and s in 0(n 2 ) so that Step III requires 
0(n 2 ) steps. Finally Step IV can be carried out in 0(n) steps so that the unranking proceedure 
requires 0{n 2 log{n)) steps. 

For an example of our unranking algorithms for C^ ± , suppose that we want to find the tree 

T G Cf ! whose rank is 600,000. In this case, (9 - 5)\S 7A = 8400 and SV i4 = 350. Thus for Step 
I, we find that 

600, 000 = 71 x 8400 + 3600 

3600 = 10 x 350 + 100. 
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Thus n = 71, r 2 = 10 and r 3 = 100. 

For Step II, first we apply the uranking proceedure for T>To,^ from section 4.2 to find 
g G T>T§£ such that rankx>T 9 5 (g) = 71. We set m' = 71. 

Then 5(1) - 1 = max{y : (|) < 71}. Note that (jj) = 56 < 71< 126 = (J). Thus 5(1) -1 = 8 
and hence 5(1) = 9. We then set m! = 71 — 56 = 15. 

Then g{2) - 1 = max{y : (|) < 15}. Note that (®) = 15. Thus 5(2) - 1 = 6 and hence 
c/(2) = 7. We then set m' = 15 - 15 = 0. 

Then g(3) - 1 = max{y : (|) < 0}. Note that Q = 0. Thus 5(3) - 1 = 2 and hence 
5 (3) = 3. We then set m' = 0. 

Then 5(2) - 1 = max{y : < 0}. Note that Q = 0. Thus 5(2) - 1 = 1 and hence 
5(2) = 2. We then set m' = 0. 

Then g(l) - 1 = max{y : («) < 0}. Note that (°) = 0. Thus #(1) - 1 = and hence 
9(1) = 1. 

Thus 5(5)) = (9,7,3,2,1). 

Next we apply the unranking proceedure for XO4 from section 4.3 to find h 6 XO4 such that 
rankxo 4 (g) = 10. We set /i(4) = and po = 10. 

Then 10 = 1 x 3! + 4 and 4 < 3! so that h(l) = 1 and Pl = 4. 

Then 4 = 2 x 2! + so that h(2) = 2 and p 2 = 0. 

Then = x l! + so that h(3) = and p 3 = 0. 

Thus (h(l), h(2), h(3), h(4)) = (1, 2, 0, 0) and h is the direction insertion sequence of the permu- 
tation cr = 2413. 

Finally we apply the unranking proceedure for TZQt^ from section 4.4 to find s = (a^, . . . , cto) 
such that rankjig 7 4 (s) = 100. 

First we set «6 = and r = 100. 

Then m 5 = a 6 = and by Table 2, £( 4 )(5,0) = 65. Thus the maximum value of s such 
that s < m 5 + 1 and s£( 4 )(5, 0) < 100 is s = 1. Thus a 5 = 1 and we set r = 100 - 65 = 35. 

Then = max({a.Q, 05}) = 1 and by Table 2, E'( 4 )(4, 1) = 55. Thus the maximum value 
of s such that s < 771,4 + 1 and sE^(4, 1) < 35 is s = 0. Thus «4 = and we set r = 35. 

Then m 3 = max({a 6 , a 5 , a 4 }) = 1 and by Table 2, £( 4 )(3, 1) = 9. Thus the maximum value of 
s such that s < 7713 + I and s£ (4) (3, 1) < 35 is s = 2. Thus ct 3 = 2 and we set r = 35 -(2 -9) = 17. 

Then m 2 = max({ae, a 3 }) = 2 and by Table 2, #'(2, 2) = 7. Thus the maximum value of 
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s such that s < m<i + 1 and sE^ (2, 2) < 17 is s = 2. Thus 02 = 2 and we set r = 17 — (2 • 7) = 3. 

Then mi = max({as, . . . ,02}) = 2 and by Table 2, .E?( 4 )(l,2) = 1. Thus the maximum value 
of s such that s < mi + 1 and sE^\l, 2) < 3 is s = 3. Thus ai = 3 and we set r = 3 — (3- 1) = 0. 

Then mo = max 

({a 6 , • • • ,«i}) = 3 and by Table 2, #( 4 )(1,3) = 1. Thus the maximum value of 
s such that s < m + 1 and s£ (4) (l, 2) < is s = 0. Thus a x = 0. 

Thus s = (0, 1, 0, 2, 2, 3, 0) which corresponds to the set partition ((1, 3, 7), (2), (4, 5), (6)). 
For Step III, we want to construct / G J-g such that 

Seq(f) = ((9, 7, 3, 2, 1), (1, 2, 0, 0), (0, 1, 0, 2, 2, 3, 0)). 

Now tt-(/) = ((1, 3, 7), (2), (4, 5), (6)) Since a = 2 4 1 3, the ordered partition vr(/) = ((2), (6), (1, 3, 7), (4, 5)) 
and the ordered partition vr + (/) = ((3), (7), (2,4,8), (5,6)). It follows that 

Table(f) = (0, 0, 0, (3), (7), (2, 4, 8), 0, (5, 6), 0). 



Thus / is specified by the following table. 



r x (i) 


r i (2) 


/ _1 (3) 


r A (4) 


r A (5) 


r A (6) 


r A (7) 


r x (s) 


r x (9) 











{3} 


{7} 


{,2,4,8} 





{5,6} 






The graph of / and the graph of ©(/) are pictured in Figure [9j 




Fi gure 9: The element of rank 600,000 in C| ^. 

The ranking and unranking algorithms for where L is a subset of [n] of size k is 

essentially the same as the ranking and unranking algorithms for t except that the decreasing 
function g E T>T n ^ is precified by L. That is, if L = {l\ < . . . < Ik}, then g = (g(l), ■ ■ ■ ,g(k)) = 
(lk,k-i, ■ ■ -,h)- Now 

l^n'l I =(n-k)\x S n - 2 ,n-k- 

Ranking Algorithm for 

Step 1 Given TeC^, construct f = G~ 1 (T). 
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Step 2. Find h G IO n -k, and s G lZG n -2, n -k such that 

5eg(/) = (( 5 (1), • ■ • g(k)), (h(l), ...,h(n-k)), (s(l), ...,s(n- 2))). 

Step 3. Use the algorithms of section 4 to compute rankxo n ^ k (h) and rankng n _ 2 n _ k (s). 

Step 4 It then follows from our analysis in section 4 that 

rank*k,L = rankxo n _ k (h) ■ S n ^ 2 ,n-k + rankTZQ n - 2t n-k(s)- (47) 

n,l 

Thus ifL = {3,5,6,8,9},n = 9, and T is the tree pictured in Figure 3, then it follows from our 
previous calculations when we computed rank^ 5 , that 

°9,1 

ran£u5,{3, 5 ,6,8,9} (T) = (7 x 350) + 46 = 2, 496 (48) 

°9,1 

Thus the T pictured in Figure 3 is the tree with rank 2,496 in (^>1 3 > 5 ' 6 > 8 ' 9 j\ 

~*k L 

The unranking algorithm for C \ again is essentially the same as the unranking algorithm 
for 0^,1 exce Pt that we do not have to find the decreasing function g since it is prespecified by L. 

Unranking Algorithm for C\ x . 

Problem: Given r such that < r < (n—k)lxS n -2, n -k, find T G \ such that rank^k,L (T) = r. 

71,1 

Step I Find n and r 2 such that 

r = n x S' n _2,n-fc + ^2 where < r 2 < S , n _ 2 ,„-fc- (49) 

Step II Find h G !O n -k and s G TZG n -2,n-k such that 

n = ranfc xc , n _ fc (/i) (50) 

r 2 = rank n g n _ 2n _ k {s) (51) 

Step III Find / G T n such that 

Seq(f) = {(g(l), g(k)), (h(l), ...,h(n-k)), (s(l), ...,s(n- 2))). (52) 

Step IV Let T = 0(f). 

For an example of our unranking algorithms for let L = {3,4,6,8,9} so that g = 

(9, 8, 6, 5, 3). Suppose that we want to find the tree T G Cg'f whose rank is 6,000. In this case, 
(9 - 5)!SV,4 = 8400 and SV i4 = 350. Thus for Step I, we find that 

6,000 = 17 x 350 + 50. (53) 
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Thus n = 17 and r 2 = 50. 

For Step II, first we apply the uranking proceedure for XO4 from section 4.3 to find h £ XO4 
such that rankxo 4 (g) = 17. We set n(4) = and po = 17. 

Then 17 = 2 x 3! + 5 and 5 < 3! so that h(l) = 2 and p 1 = 5. 

Then 5 = 2 x 2! + 1 so that h(2) = 2 and p 2 = 1. 

Then 1 = 1 x l! + so that /i(3) = 1 and p 3 = 0. 

Thus (h(l), h(2), h(3), h(4)) = (2, 2, 1, 0) and h is the direction insertion sequence of the permu- 
tation cr = 3 4 2 1. 

Then we apply the unranking proceedure for TZQt^ from section 4.4 to find s = (as, . . . , ao) 
such that rank-Rg 7 4 (s) = 50. 

First we set «6 = and r = 50. 

Then 777.5 = «6 = and by Table 2, £^(5,0) = 65. Thus the maximum value of s such 
that s < 7775 + 1 and sE^(5, 0) < 50 is s = 0. Thus 05 = and we set r = 50. 

Then 7774 = max{{a 6 ,a 5 }) = and by Table 2, £( 4 )(4,0) = 10. Thus the maximum value 
of s such that s < 7774 + I and s£ , ^(4, 1) < 50 is s = 1. Thus a 4 = 1 and we set r = 50- 10 = 40. 

Then 7773 = max({aQ, a^, 04}) = 1 and by Table 2, E^(3, 1) = 9. Thus the maximum value of 
s such that s < 7773 + 1 and s£( 4 )(3, 1) < 40 is s = 2. Thus a 3 = 2 and we set r = 50- (2-9) = 22. 

Then tt7 2 = max({a 6 , a 3 }) = 2 and by Table 2, £( 4 )(2, 2) = 7. Thus the maximum value of 
s such that s < m 2 + l and s£( 4 )(2,2) < 22 is s = 3. Thus q 2 = 3 and we set r = 22 - (3 • 7) = 1. 

Then mi = max ({a 6 , . . . , a 2 }) = 3 and by Table 2, £( 4 )(1, 3) = 4. Thus the maximum value of 
s such that s < mi + 1 and si?( 4 )(l, 3) < 1 is s = 0. Thus a± = and we set r = 1. 

Then m = max({a 6 , . . . , ai}) = 3 and by Table 2, £^(0, 3) = 1. Thus the maximum value of 
s such that s < m + 1 and s.E( 4 )(l, 2) < 1 is s = 1. Thus a\ = 1. 

Thus in our case, s = (0, 0, 1, 2, 3, 0, 1) which corresponds to the set partition ((1, 2, 6), (3, 7), (4), (5)). 
For Step III, we want to construct / G Jg such that 

Seq(f) = ((9, 8, 6, 5, 3), (2, 2, 10), (0, 0, 1, 2, 3, 0, 1)). 

Now -K-(f) = ((1,2, 6), (3, 7), (4), (5)). Since a = 3 4 2 1, the ordered set partition tt(/) = 
((4), (5), (3, 7), (1, 2, 6)) and the ordered set partition vr+(/) = ((5), (6), (4, 8), (2, 3, 7)). It follows 
that 

Table(f) = ((5), (6), 0, (4, 8), 0, 0, (2, 3, 7), 0, 0). 
Thus / is specified by the following table. 
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r x 0) 


r i (2) 


r A (3) 


r A (4) 


r A (5) 


/ _1 (6) 




/ _1 (8) 


/ _1 (9) 


{5} 


{6} 





{4,8} 








{2,3,7} 









The graph of / and the graph of 0(/) are pictured in Figure 10 

" V i 



1 = 




3 8 



0(f)-- 




Fi gure 10: The element of rank 6,000 in Cg'| \ 



Since the ranking and unranking algorithms for C n \ are essentially a subset of the ranking 
and unranking algorithms for 1 , it follows that both the ranking and uranking algorithms for 



C n \ require 0(n 2 log(n) steps. 
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